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The bivariate probability distribution for two composite meteor signals displaced in 
time is derived theoretically using the Markoff statistical combination technique. Both 
the effects of numerous, small meteors and the residual reflections from infrequent large 
meteors are treated simultaneously. For the case of exponential decay of component signal 
spikes which are themselves distributed as the inverse square of their initial amplitudes, we 
find that the joint probability that a composite signal ki is observed at time ti and R2 at 
t'z, seconds later, is given exactly by the following expression involving elliptic functions: 
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where a = Qvr} and v is the average rate of occurrence of meteor signal spikes of all sizes and 
7f is the exponential decay time of each initial spike. This result reduces to the usual limit- 
ing forms in the case of r very large or very small relative to the decay time -n. 

1. Introduction 

The reflection of VHF radiowave signals from meteor trails in the A-region of tlie ionosphere 
is interestin<>: for several reasons. It occurs both as radar hackscallerintj^ from the trails and 
as propagation over ol)li([U(» paths ranging out to 1,500 km. The latter is significant because 
of its communication opportunities, and in tliat context it is important to know as much about 
the signal structure as possible. The composite meteor signal is the result of reflections from 
both large and small meteor trails, all in various stages of decay. The large meteors are easily 
recognized as individual signal spikes in amplitude versus time records. However, the far 
greater iiuml)er of small, indistinguishable meteors also make an important contribution to the 
composite^ signal. The basic problem of unch^rstanding the signal structinv is to treat the con- 
tinuum of trail sizes simultaneously. 

Th(^ ])robability distribution for the envelope of a meteor signal composed of reflection 
contributions from many trails of various sizes was derived theoretically in an earlier paper 
[Wheelon, I960]. Both the effects of numerous, small meteors and the residual reflections 
from infrequent, large meteors were treated simultaneously. For the particular case of expo- 
nential decay of initial spikes which are themselves distributed as the inverse square of tlieii' 
amplitudes, we found that the probability that the composite signal amphtude sbould c^xcccmI 
a prescribed level r is given by: 

P(R>r)=^ \^^,> (1.1) 

where v is the average rate at which echoes of all sizes occur, 77 is the charact(M;istic decay time 
of each meteor spike and is set by the diffusive decay of the trail itself. Expression (1.1) be- 
haves like a Rayleigh distribution for small amplitude margins r. For the larger, less likely 
signals r, it agrees with the result predicted by elementary analysis of isolated meteor reflec- 
tions in various average states of decay. 
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Essentially the same result was established by Bain [1960] in Britain and published almost 
simultaneously. His treatment goes to somewhat greater lengths to remove a potential diver- 
gence for small signal values by renormalizing the constant of proportionality Q in the assumed 
distribution {Qjp'^) for the strength of individual spikes p in terms of an effective amplitude 
cutoff e. However, it would appear that this procedure is not required, since both (1.1) and 
the amphtude distribution corresponding thereto, 
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are well behaved for small signal amplitudes. It seems that the statistical combination of 
individual signals discriminates against the small signals, and that the weighted sum is insensi- 
tive to the error made in trying to extend the function Q/p^ (which was curve-fitted to the 
experimental data) to smaller values of p. 

In point of fact, it is the behavior of this assumed distribution for large p which causes 
trouble, in the sense that (1.2) does not possess finite moments of any order. This is due to 
its behavior for exceptionally large signals, 
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which can, in turn, be traced to the mild behavior of Q/p'^ for large p. 

The purpose of the present paper is to extend our knowledge of meteor signal structures 
by deriving the bivariate probability distribution for composite meteor signals displaced in 
time by an interval r. Such a result is evidently important in communication applications, 
since the time that a meteor scatter circuit is open is related to the interval during which all 
signals are above a specified threshold. The correlation of signal amplitudes between two 
instants is another measure of communication capacity. All of these basic features of a meteor 
circuit are derivable from the joint probability, 
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that the composite signal assumes a precise value 7?i at time t, and becomes exactly i?2 at a 
time t2, T seconds later. The probability density W {R1R2T) will be derived in this paper. A 
third paper in this series will use this result to calculate communication characteristics of 
meteor scatter circuits. 

2. General Expression for Bivariate Distribution 

As in the first paper of this series, we shall use the Markoff method [Chandrasekhar, 1943] 
to calculate the probability distribution for composite meteor signals. First, a word about 
notation and convention. Consider the typical signal history shown in figure 1, where time is 
run positively toward the past for analytical convenience. The two instants at which we wish 
to estimate the signal probability distribution are denoted by ti and #2; although ^2 can and 
will later be chosen as the time origin. A large but finite interval T is chosen in which A^ 
meteor bursts are assumed to have occurred. A/' is a statistical quantity, whose mean value is 
vT, where v denotes the average rate of meteor occurrence. The limit as T goes to infinity 
will be taken later in the analysis, after convergence of certain integrals is assured. 




Figure 1. Typical succession of meteor spikes with 
initial amphtude p, indicating measuring instants 
ti and t2, separation time r, and total interval T. 
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The basic statistical problem is to calculate the probability density for observing com- 
posite signals /?i and /?2 at U and ti respectively. These total signals are compounded from 
the residual signal elements of all prior meteor spikes in the interval T. All signals received 
prior to to contribute in some residual way to both Bi and i?2. On the other hand, spikes 
received after ^2 but prior to ti contribute only to Ri. Of course, the individual spikes con- 
tribute to i?i and B2 in different measure, because of differential amplitude decay and phase 
relationships. However, one can write explicit expressions for the two total signals, 

Ri=JlvMti-U), (2.1) 

in terms of the discontinuous function 

, KO, (2.2) 

which has been introduced for analytical convenience. This function allows one to write 

both i?i and Ri as formal sununations over all A^ initial spikes pi, yet ensures that those spikes 

received after the individual measui'ing events arc not actually included. This is essentially 

a bookkeeping device, but is of considerable utility in organizing the subsequent analysis. 

-> -» 

The joint probability distribution for Ri and R2 is given by Markoff's general metliod 

[Chandrasekhar, 1943] as a double Fourier transform of the characteristic function Bijci.ko) 

Wiliji^ = (27r) -' Ll'k, P/2^26'^^-^i+'2-«2) B {k,h) (2.3) 
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where 



The brackets indicate that one is to average over the three statistical features of the individual 

signals: (1) the distribution of initial pulse heights j^t and their random phase, (2) the prob- 

-> 
ability of finding a spike jh at the instant tt, and (3) the actual number of echoes N in the 

fixed interval T. Because the individual meteors are statistically independent of one another 

(i.e., i), one can write 

B{k^k2)^((exp i[k^-^\(t-t,)+k2'p\(t-t2)]YS . (2.5) 
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Since the train of meteor echoes apparently forms a Markoff process (no sense of history) of 
small probability, one can argue that the probability of observing exactly N spikes in the 
interval T should follow a Poisson distribution. 

PiN/T)^^ e-'^. (2.6) 

The average over N of (2.5) is computed thus: 

B{kik2)=J2 ^'^^ g-.T /^^fci.px(^-.,)-fa2./^xu-^2M _^ 
N=o iV! \ /t,p 
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One can write the remaining averages over p and t in terms of the probabihty yip,t) that a 
single meteor echo occurs at time t and produces a vector signal p in the receiver. 

lCj^-p\(t-t{) + ik2-p\(t-t2) \ -- 1^2^ I ^^yfj)-f^\^iki.p\it-t^) + ik2.p\it-t2) (2 8) 

/t,p J Jo 

Since the individual meteor echoes can occur with equal probability anywhere in the inter- 
val T, 

-> 1 -> 

The initial echo spikes p are randomly phased, since the distance from the transmitter to the 
individual meteor trails is a completely random variable. Hence, 
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where 6 is an arbitrary phase reference, which we take to be the angle between p and ki. D(p) 

-> -^ 

represents the distribution of initial pulse heights. If <^ denotes the angle between ki and k2^ 

one can use the equivalence (2.8) and the above to rewrite expression (2.7) for B(kik2) as 

follows: 

B{kik2)=exp-^ f dpD{p) {^ dd{ ^^[i_g':^i/'cos0x(^-/i)giVeo8(0+<^)x(^-^2^]^ (2.9) 
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The angular integration can be done using the standard formula: 

P" dde'^'' ^^° ^+'' ^"^ ^^ = 27rJo(V^'+t') (2.10) 

giving: 
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where we have taken the limit of T going to infinity, since the difference quant it}- in 
braces is now finite at the upper limit. 

To make further progress, one must divide the time interval and introduce explicit expres- 
sion (2.2) for the discontinuous time functions X (^— ^i) and X (^—4) • At this point in the calcula- 
tion, it is convenient to choose ^i = and^2 = 'r, sinceonly the time difference (^i— ^2) = ^ is relevant 
to a stationary problem. Dividing the time interval into two segments: 0</<r and r<^<oo, 
allows one to write equation (2.11) out as follows: 

S(ti,/:2,0) = exp-.r^/:pD(p)/ Jjr//[l-Jo(^2P6-^/^)] 

+ I dt[\-J,{pe-"^^|kl+2k2k, cos ct>e^i^+k\e'^'^)] V- (2.12) 

3. Distribution for Inverse Square Law 

The remaining integrals can be done in closed form if the distribution of initial pulse heights 
is assumed to have the form, 

Z>(2,)=|. (3.1) 
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This expression is jiiore than an analytical convenience, in that experimental data on the 
distribution of individual pulse heights is well fitted by eq (2.13). The major assumption 
involved in using (2.13) is that the same law extends down to the small meteors which cannot 
be distinguished as individual spikes. This assumed distribution of initial pulse heights is not 
reall}^ an acceptable form in that the normalization integral 
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diverges at the lower limit, rather than approaching unity. Strictly speaking, this implies that 
the average rate of meteor signal occurrence v is infinite, which, in turn, invalidates the assump- 
tion of the Poisson distribution (2.6). The cutoff procedure introduced by Bain [1960] remedies 
this deficiency in a formal way. However, the divergence of this law for small p does not affect 
the final composite signal, as noted earlier. Its relatively slow decrease for the very large, 
exceptional rare meteors causes the real trouble. 

The double integrations remaining in the evaluation of this characteristic function B 
{kik-i) in eq (2.12) can be done by interchanging the order of p and / iiitogi'at ion mid setting 

and 



v=pe-'^^^Jkl+2k2kl cos </>6^/''+/:2g2r/, (3^2) 

in the first and second integral, respectively. This gives for (2.12) the follow ing: 

+VAri+2A-2^i cos <t>e^'^+k\e'^'^' C dte-"^ T J [1-JnC^Ol }• (3.3) 
The integrations are now uncoupled and can be done by noting that 

r"$[i-jo(i-)]=i. (3.4) 
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The final expression for the characteristic function becomes. 



B{kuk2,(t>)=^'^V-^vQ [k2{l-e-''^)+-yJkle-^'''^+2k2k, cos 4>e-'^''^+k~\' (3.5) 

The bivariate probability density is calculated from this expression as the double Fourier 
transform of eq (2.3) 

W(i?J?2)=-Jy4 \d'k, { d'k2e''^'y'^^'^'^'^-^^^ e-'^^2^'-'~'^'\-^^ (3.6) 

where for notational convenience we have now set 

c=vy)Q. (3.7) 

The integrals in eq (3.6) can be performed most readily if one makes the following linear 

vector transformation: 

-> ^ -^ 

so that 

The / and ki integrals are now separated and can be done analytically by using the intra-vector 

angular definitions exhibited in figure 2. With these conventions, one can write out the four- 

— * -* 

fold integration above as angular definitions for signal (/I'l) and transform {ki) vectors, all 
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Figure 2. Intra-vector angular definitions. 



with respect to an arbitrary phase reference vector /. 
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Using the integral result (2.10), we find the final expression for the bivariate vector probability 
distribution to be 
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One can verify that this expression is normalized to unity by integrating over the compo- 
nents of Ri and i?2. 
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The distribution of signal amplitudes is obtained from (3.8) by integrating over the phase 
angles ^ and 6, and leads to the following expression involving complete elliptic functions. 
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The expected limiting forms of this result emerge if one assigns special values to the time dis- 
placement T. For example, if the time difference between the two measuring instants is large 
compared to the decay time of the individual signal specified n, one has 
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This is just the product of the individual distributions (see eq 1.2) for i?i and R2 derived in 
the first paper of tliis series [Wheelon, I960]. The opposite extreme is somewhat more subtle, 
but can be extracted by taking the limit as r goes to zero. 
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Using the following limit definition of the Dirac Delta function, 

3(x)=-lim( ^-f— 2 ) 
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we have 
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This indicates that as the time displacement goes to zero, the two signals should coalesce as 
prescribed by the Delta function and their distribution be described by the previous result 
for the single time expression. Note that the factor Ki in the Jacobian for polar coordinates 
is just cancelled by the (hMioniinjitor tci'iii. 
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